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1. Introduction 

The calculations of multi-particle processes at next-to-leading order (NLO) level are 
particularly important for Large Hadron Collider (LHC) and International Linear 
Collider (ILC) physics, especially in hunting for Higgs boson(s) and New Physics 
(NP) signatures. Searching for the Higgs boson is one of the central tasks required 
for the operation of LHC. In order to solve some difficulties encountered in the Stan- 
dard Model (SM), many New Physics models have been introduced. Most notable 
of them are the various supersymmetric (SUSY) models. The Minimal Supersym- 
metric extension of the Standard Model (MSSM) [J], g, |, |, § §, 0, §, the simplest 
one, provides many new physics benchmark points to experimentalists to search the 
SUSY signatures at the LHC. However, both hunting for the Higgs boson(s) and 
searching for the New Physics signatures are often limited by overwhelming back- 
grounds. Cascade decays from heavy SUSY particles often result in multi-particle 
final states. Therefore, precise theoretical results of multi-leg processes are urgently 
needed with the enhancement of accumulated data at the LHC. To accomplish this 
task, a significant number of multileg processes should be calculated up to the NLO 
accuracy in SM and NP. Nevertheless, this task can be very challenging. 
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In general, a NLO calculation includes real-radiation corrections and virtual 
corrections with renormalization. The real corrections, without any loop integrals 
at NLO level, can be accomplished with many efficient algorithms, such as Berends- 
Giele recursion relations^, Britto-Cachazo-Feng-Witten (BCFW) recursion relations 
[0, 0, etc. Programs such as MADGRAPH0, |IJ, COMPHEPjlJ, AMEGIC++[|15j 
, ALPHA |T§, [T7|, [Tg|, HELAC[0, [2T| and COMIX[gT[] are all efficient tree-level matrix 
elements generators. Accompanied with phase space slicing method P2| or subtrac- 

|0, EyJ, the automation of real-radiation part 



tion terms [|23|, |2J, |25|, |2§ £7|, |28 
seems much straightforward to be realized [|32|, ^7 , p3| 0. On the other hand, the 
automatic calculation of one-loop integrals is also a feasible task nowadays. Actu- 
ally, the difficulties in one-loop calculation are relevant to the simplification of the 
lengthy Dirac structures and the reduction of one-loop integrals to standard master 
integrals. The latter issue can be resolved by many existing one-loop integrals re- 
duction algorithms, such as Passarino-Veltman reduction procedure [B5L |36|, B7|] . The 



former one encounters difficulties in making the expression more compact, because 
one should treat the Dirac matrices in d = 4 — 2e dimensions when using dimensional 
regularization and dimensional reduction. 

Recently, we have witnessed a new evolution of NLO techniques. Automatic 
one-loop calculations have become a feasible approach after several new and efficient 
algorithms are proposed. Some of the most notable methods are the Unitarity [[38, 
|39| , f4C| |4l| , [42], [|3| based techniques such as the Ossola, Papadopoulos, and Pittau 
(OPP) reduction method P4L HR [46], E7|, [48], |49|], by reducing the computation of one- 



loop amplitudes to a problem of a tree level calculation. This method are able to 
control the computational complexity efficiently with remaining tree level recursion 
equations. However, when applying in 4 dimensions in OPP, one can extract the Cut 
Constructible(CC) part of the amplitudes, while the left piece rational terms should 
be derived separately [[SOR . With the known rational part R, efficient simplifications of 
the lengthy Dirac structures in 4 dimensions are possible at the Feynman amplitude 
level, which also make it possible to do a complicated one-loop computation about 
multi-particle processes in other reduction algorithms [plj, [ 



53 



Fortunately, the rational part R (in OPP approach, it is called R2) is proven to 
be guaranteed by the ultraviolet (UV) nature of one-loop amplitudes J55], 



,i.e. the 

only origin of R we considered here is from a combination of 0(e) part in numerator 
of a loop integral and its UV divergence term 0{-^) . Since the UV poles, in contrast 
to infrared divergence ones, do not depend on kinematical properties of external legs 
such as on-shell relations, one can establish the Feynman rules respect to the effective 
vertices. This fact also ensures four external legs enough. 

The complete Feynman rules for R in SM under anticommutating 75 strategy 
in the 't Hooft-Feynman gauge, i?g gauge and Unitary gauge have been calculated 
[ j5"5| , |56| , ISTf . Besides, their package R2SM written in FORM is also available [ |58|j . 
The Feynman rules for the R in SM under the 't Hooft-Veltman 75 scheme have 
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been calculated by us |p9| . Moreover, some simplifications to extract rational terms 
were suggested recently |Kj ^T] . For supersymmetric amplitudes, one can track to the 
Cachazo-Svrcek-Witten Fenyman rules p2|, |63| , |64]| to make some simplifications in the 



calculation of gluonic amplitudes. In the present paper, the complete set of Feynman 
Rules for rational part R of QCD corrections in the MSSM are calculated at the one- 
loop level with two 75 strategies. All of these Feynman rules can be implemented in 
NLO matrix element generators like MADLOOP^, and HELAC-NLO[[T7| and 
also be useful in the development of FEYNRULES |68H or in any other methods such 



as Open Loops [§9| or GOSAM [0, |7| 



The organization of the paper is as follows. In Section 2, the origin of rational 
part R is recalled. The dimensional regularization schemes and 75 schemes are fixed 
in Section 3. The complete set of Feynman rules are listed in Section 4. Finally, we 
make a conclusion. 



2. Origin of Rational Part 

In dimensional regularization procedure, one should treat integral momentum q in 
d = 4 — 2e dimensions to maintain the gauge invariance. A generic A-point one-loop 
(sub-) amplitude reads as 

where 

D k = (q+Pk) 2 - {m k ) 2 , q = q + q. (2.2) 

Here a bar denotes d dimensional objects while a tilde means something living in 
d — 4 dimensions. 

The rational part R comes from the division of the above numerator N(q) into 
a 4-dimensional part and a (d — 4)-dimensional part 

N(q) = N(q) + N(q 2 ,e iq ), (2.3) 

Apart from the process dependent N(q), the remaining part N(q 2 , e, q) = N(q) — N(q) 
gives rise to rational part R 

R = [ dtqMl^L, (2.4) 

To clarify the division and to avoid the possible ambiguities, we split d = 4 — 2e 
dimensional objects in the tree like Feynman rules as follows 

?A = fy + ?£, 7a = Tm + Ta 3 9js» = 9 V v+9jis>- ( 2 - 5 ) 
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The effective vertices for R can be obtained from all possible one-particle irreducible 
Green functions, which is enough up to 4 external legs in the MSSM QCD, which 
is a renormalizable model. It's also worth reminding that the rational part R is 
not gauge invariance independently. Actually, in OPP framework, Ri + R 2 is gauge 
invariant and preserves all the Ward identities of the theory [|56|j. 



3. Dimensional Regular izat ion Schemes and 75 Schemes 



In this section, the proposed regularization schemes and 75 schemes adopted in this 
paper are reviewed and an example of calculating rational term R is given. 

First, the Feynman rules in two dimensional regularization schemes that main- 
tain the advantages of the helicity method for loop calculation are reviewed. These 
schemes require all the quantities of external legs to be in four dimensions. The 
supersymmetry-preserving Dimensional Reduction (DRED)[^j has been proven to 
be equivalent to Dimensional Regularization (DREG) ff3j . Thus, the Four- dimensional 
Helicity (FDH) f7|, [75|, [?7| and 't Hooft-Veltman (HV) [78] schemes were chosen 
in our results. In the FDH scheme, the only d dimensional object is the integral 
momentum q, while in the HV scheme all the internal (unobserved) quantities such 
as the polarization of internal vectors are d dimensional objects, i.e. 



R 



HV 



d,N(q 2 ,e,q) 



d d q 



D X D 2 ...D 



N 



R 



FDH 



d d q 



_N(q 2 ,e = 0,q) 



D 1 D 2 ...D 



(3.1) 



N 



When using dimensional regularization in both schemes, the common object to be 
analytically continued from 4 dimensions to d = 4 — 2e dimensions is 



q 2 ->• q 2 + q 2 . 



(3.2) 



Because of the nature of external legs in four dimensions, the contribution of q to R 
has to be in q 2 forms. 

As we known, there are two famous algebraic self-consistent 75 regularization 
schemes existing in the market. The first class, called 't Hooft-Veltman (HV) 75- 
regularization scheme, was proposed by 't Hooft and Veltman and systematized by 

It is the first proven 



Breitenlohner and Maison TO gO, Bl], g2, |3, 8J, g5, 



consistent 75 scheme [|82| ^3| ^4[ £5|, It favors a non- vanishing anti-commutator 
{ihilfi} 7^ 0. This scheme distinguishes 4-dimensional and (d — 4)-dimensional ob- 
jects to create correct spurious anomalies. Due to the non-anticommutation relation, 
the tree-like Feynman rules in the theory were modified as symmetric forms []87| in 
our calculations. The other practicable scheme was introduced by Kreimer, Korner, 
and Schilcher (KKS), et al. PSI EO, BD1. Instead of violating the anticommutation 
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equation, they chose to preserve the anticommutation relationship {75, 7m) = and 
prevented cyclic property in Dirac traces to avoid algebraic inconsistency. A projec- 
tion on four-dimensional subspaces is needed in their redefinition of trace operation: 

^(757^-7^) = MP(75)7w-"7fc,J, ( 3 - 3 ) 

where V{^) = \ x e ^2^4 7 MI l m ! 7 /i3 7 At4 with Lorentz indexes of £ WM2jU3At4 all in 4 
dimensions. To obtain the correct result in this scheme, a "special" vertex called 
"reading point" should be identified in the loops and all the 75 in Dirac structures 
should be removed to the vertex before performing projection. Generally, different 
treatments of "reading point" in each Feynman graph may generate a discrepancy 
proportional to the totally antisymmetric Levi-Civita tensor e. In some computa- 
tions of specific processes, different intermediate results may be derived if different 



regularization schemes are selected [91], |9l|. However, physical results should be the 



same in different regularization schemes as long as the unitarity is kept in the theory 



after some necessary renormalization |p3|, |94], |95], £)J| [57], pq| . 

To describe the procedure more clearly, an illustrative example of deriving R 
from the gluon self-energy in MSSM QCD is provided. The corresponding Feynman 
diagrams are shown in Fig.|l]. The contribution of diagram (a) in Fig.|I] is vanishing 
because it is a massless tadpole. For (b) — (d), the internal scalar loops cannot give 
a non-vanishing contribution to the R of gluon self-energy because the vertices are 
always contracted to external polarization vectors. The nontrivial contributions are 
only from the last three diagrams. For the quark loop with two external gluons, the 
numerator can be read as 

N(q) = -7^Y Tr[r(l + mQ)r(l + ^ + mQ)] ' (3 ' 4) 

where external indices \x and v have been taken in four dimensions. After performing 
some Dirac algebra, one arrive at 

N(q 2 ) = IjiyJf^Q 2 - (3.5) 

Integrating it with the help of any one-loop integral reduction algorithms, this quark 
loop contribution can be obtained. The last two diagrams are gluino loop and gluon 
loop. The same basic procedure is employed to deal with these loops. 



4. Feynman Rules for the Rational Part 

In this section, the complete list of the effective MSSM QCD vertices contributing 
to R in the 't Hooft-Feynman gauge are given. We use FEYNARTS [f59[] to gener- 
ate all of the Feynman amplitudes. Therefore, the Feynman rules are following the 
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Gl ^mSmfmrnu^ g\ gi ^mrnm*mum^ g\ 

0) (b) 



4f&6T^ Gl 



(c) 



Gl ^rrj4 



(d) 




conventions of Refs.pM 101, 102, 103 



In particular, fermion chains are concate- 
nated following the algorithm in Ref. [|104|| , which works also for Majorana fermions 
and the fermion- number- violating couplings. Two parameters Xhv an d g5s are in- 
troduced in the formulae to denote the different dimensional regularization schemes 
and 75 schemes as in our previous paper [[59|. Here, Xhv = 1 corresponds to the 't 
Hooft-Veltman regularization scheme and Xhv — to the Four Dimensional Helicity 
scheme, while g5s = 1 corresponds to the KKS 75 scheme and g5s = — 1 to the HV 



75 scheme. Our notations are as follows: 

Li = e, L 2 = n, L 3 = r, 
L 4 = v e , L 5 = z/ M , L 6 = v T , 
Qi = u, Q 2 = d, Q 3 = s, 
Qa = c, Q 5 = b, Q 6 = t, 
L\ : i = §1, Lip = e 2 , L>2 y i = jjLi, 

L*2,2 = fj>2, ^3,1 — T\, L 3j2 = T~2, 

L 4 = D e , L 5 = v^, Lq = u T , 

Ql,l = Ml, Ql,2 = «2, ^2,1 = dl, 
£?2,2 = d 2 , $3,1 = Si, (^3,2 = S2, 
£?4,1 = Ci, (^4,2 = C 2 , (55,1 = ^1, 

Q5,2 = b 2 , Q&,l=t u Q6,2=h- (4.1) 

In addition, 

ei = e, e 2 = u, e 3 = r, 

^1 = f e , ^2 = ^, ^3 = K,-, 

Ui = u, U 2 = c, C/ 3 = t, 
D l = d,D 2 = s, D 3 = b, 
h,i = ei, e 1>2 = e 2 , e 2j i = //i, 
e 2 ,2 = £2, e 3) i = fi, e 3i2 = f 2 , 

^1 = Z>e, ^2 = ^, ^3 = S,-, 
^1,1=^1, U lt2 =U 2 , U 2 ,1=C!, 

U 2:2 = c 2 , ?7 3i i = ii, C/ 3)2 = t 2 , 
D 1:1 = d 1 , D h2 = d 2 , £>2,i = h, 

D 2 , 2 = s 2 , ^3,1 = &i, ^3,2 = 62. (4.2) 

iV c denotes the number of colors and are CKM matrix elements. In MSSM, 

the following rotation matrices should be introduced too 



Qt,i \ _ R ( Qi,l Wcr,,i \= RU ( Ui, L 



Ql,2 J \Ql,R J \Ul,2 J \Ul,R 

Dl,l \ „n ( Dl,L \ ( h,l \ „ r ( h,L 



, RDi ~ h±J , =Mi, 
A,2 / \ D l,R J \ e l,2 J \ e l,R 

Besides, matrices RN, CR,CL are the neutralino mixing matrix and the right-handed 
, left-handed chargino mixing matrices respectively, which is also following the con- 
vention of FEYNARTS. We denote the element of matrix M as M^, M id or M (ij) . 
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To make the result readable, we introduce some notations for the following frequently 
used summations 



2 

SRl(l,m) s i,s2 = / ] Rj(sl,k)Rm,(s2,k) = [RmRl) , 
'—^ \ / s2,sl 

k=l 

SR2(l,m) s i )S 2 = R*i^ s i^)Rm,{s2,2) + Rl,(sl,2)Rm,(s2,l), 
2 

SRUl(l, m) sl , s2 = V RU; [sl k) RU m ^ k) = (RU m RUl) , 

k=l 

SRU2(l,m) s i :S2 = RUi^ sl ^RU m ^2) + ^^*( s i,2)-^m,(s2,i)j 
2 

m) sl ,,2 = V RDl (shk) RD mXsm = (RD m RDj) , 

— V / s2,sl 

fc=l 

SRD2(l,m) sltS2 = RDi^ sl ^RD m ^ S 22) + RD*,(si,2)RD m ,(s2,i)- 

Q~ ee i^ 5 - and f2 + = -4p£. are two chiral projection operators with g s ,e coupling 
constants in QCD and QED respectively. For brevity, tp = tan/3,c^ = cos (3, sp = 
sin /3, c a = cos a, s a = sina,^ = cos 9 w ,s w = sin9 w . Parameters /i and Aq 1 are the 
Higgs-doublet mixing parameter and the soft breaking A-parameters. In some cases, 
we also use notation A = B(a — > (3), which means the expression A are obtained 
from B by all of the a in B replaced by the f3. 

4.1 Effective Vertices in Pure MSSM QCD 

In this section, we give the complete list of the non-vanishing R effective vertices in 
pure MSSM QCD, with all external and internal legs MSSM QCD particles. 

4.1.1 Pure MSSM QCD effective vertices with 2 external legs 

All possible non-vanishing 2-point vertices in pure MSSM QCD are shown in FigfJ. 



Squark-Squark vertices 



There is no generation mixing in these vertices. The corresponding effective vertex 
is shown in Fig.|2| (a) with their expressions as follows 



Vert(gj iSl ,Qi, s2 ) 



24tt 2 
1 — g5s 



C F 5 ij 8 lm {(3ml + 3m 2 Qi -p 2 ) 

s l,s2 



1 + g5s 



SR2(l, l) s i t s2 



SRl(l, /)sl,s2 

(4.3) 
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= VeTt(Q{ sl ,Q m , 2 ) 



(a) 



Gl ^~0~07PQWC1T Gl = Vert(G^) 



(b) 



O 



= Vert(Qf,Qj 



(c) 



G a ""m^QFmm^ G b = Vert(G",G 6 



(d) 

Figure 2: All possible non-vanishing 2-point vertices in pure MSSM QCD. 



Gluon-Gluon vertex 



The corresponding diagram is shown in Fig.|2| (b) with its expression as follows 



Vert(Gj;, Gl) = C A 5°» 



Q 



(4.4) 



Quark-Quark vertices 



The corresponding diagram is shown in Fig.|2| (c) with its expression as follows 



Vert(Qj,Q 



m) 



16tt 2 



C F S ij Si m (-p' + 2m Ql )\ HV . 



(4.5) 
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= Vert(G a ,Oi,^„,J 



= Vert(G;,QI iSl ,0i, s2 ) 



= Vert(G° G b , G c ) 



Figure 3: All possible non-vanishing 3-point vertices in pure MSSM QCD. 

Gluino-Gluino vertex 

The corresponding diagram is shown in Fig|| (d) with its expression as follows 

Vert(G a , &) = ^C A 6 ab (-jj + 2m d ) \ HV . (4.6) 

4.1.2 Pure MSSM QCD effective vertices with 3 external legs 

All possible non-vanishing 3-point vertices in pure MSSM QCD are shown in Fig|| 



Gluon- Quark- Quark vertices 

The corresponding diagram is shown in Fig|| (a) with its expression as follows 



Vert(G^, Q\, Ql) = ^ T£ 5 lm (CJT + 7m , (4.7) 
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where 



C_ = (1 + Xhv) Cf + y ( R i R > 
C + = (1 + \ hv )C f + ^-(r\R^^ 



(4- 



Gluon-Gluon-Gluon vertex 

The corresponding diagram is shown in Fig.|3] (b) with its expression with only Nf 
flavor quarks as follows 



Vert(G*, Gl, G c p ) 



g 3 s N c (15 , , , 2Nf i fabc 



48vr 2 V 4 



— + A w + 



Nr. 



f abc V^ p ( Pl ,p 2 ,p 3 ), (4.9) 



where 



Vp(pi,P2,Ps) = g^{p2-Pi)p + 9up(jPs-P2)^ + 9p^(jH-Pa)u' (4.10) 

Gluino-Quark-Squark vertices 

The corresponding diagram are shown in Fig.| (c, d) with their expressions as follows 



Vert(G a , Ql = % Tj* <W (CLfl" + , 



7T 



Vert(G a , Q{, ) = 7™ ^ m (C^" + CT + ) , 



7T 



(4-11) 



where 



32 v / 2iV c 



J*,(.,i) - 4(1 + X H v)Ri,(s,2) + 4(1 + A HV )i2, 



-R 



1,(8,1) 



c 4 



32v^iV c 



-R, 



l,(s,2) 



1 — g>5s 
2 

Rl,{s,2) 
1 — g5s 



2./ l , / . ,.„.,) ( R}Ri) 2i + + 2 2 Ri,{ S ,2) (rIr^ | 



1 — #5s 



4(1 + Affv)i^ 



1 + g5s 



4(1 + X HV )R, 



l,(s,2) 



N 



2R 



'•(■-■2) l- R l-^) 12 + + 2 1,1 



■i,{s,i) i -Rz -R; 



1.2 



(4-12) 
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Gluon-Squark-Squark vertices 

The corresponding diagram is shown in FigFJ (e) with its expression as follows 



Vert(GlQl 1 ,Ql yS2 ) = ^T«5 lm 



6iV c 2 - 1 

~lmT dsl ' s2 



(SRl(l, l) sl , s2 - SR2(l, l) sliS2 )\ ( Pl - pa) M . (4.13) 



Gluon-Gluino-Gluino vertex 

The corresponding diagram is shown in Fig.|| (/) with its expression as follows 



Vert(G« ; G\ G c ) = ^ (C-fl" + C + Q + ) 7m , (4.14) 

where 



= 2N C (Tr(T a T b T c ) - Tr(T a T c T b )) (1 + \ HV ) 

Tr(T a T b T c ) J2( R l R i) -Tr(T a T c T b ) ^ (#J#, 
\ l ' i 

C + = 2N C (Tr(T a T b T c ) - Tr{T a T c T b )) (1 + Xhv) 



2.2 



+ ^Tr(T a T b T c ) ^(R\R^^-Tr(T a T c T h ) ^ (4-15) 



4.1.3 Pure MSSM QCD effective vertices with 4 external legs 

All possible non-vanishing 4-point vertices in pure MSSM QCD are shown in Fig^. 



Gluon-Gluon-Gluon-Gluon vertex 



The corresponding diagram is shown in FigfJ (a) with its expression with only Nf 
flavors quarks as follows 




(c) 

Figure 4: All possible non-vanishing 4-point vertices in pure MSSM QCD. 



Vert {G% G b u , G% G d a ) = ^ (d g^g pa + C 2 g w g va + C 3 g^g vp ) , (4.16) 

where 



d = Tr({T a , T b }{T c , T d }) (IIN C + 2X HV N C + 6N f ) 

- (Tr(T a T c T b T d ) + Tr(T a T d T b T c )) (22N C + AX HV N C + 10N f ) , 
C 2 = C 1 (b++c) = Tr{{T a , T c }{T b , T d }) (UN e + 2\ HV N C + 6N f ) 

- (Tr(T a T b T c T d ) + Tr(T a T d T c T b )) (22N C + A\ HV N C + WN f ) , 
C 3 = d(6 d) = Tr({T a , T d }{T c , T b }) (UN C + 2X HV N C + QN f ) 

- (Tr(T a T c T d T b ) + Tr(T a T b T d T c )) (22N C + A\ HV N C + 10N f ) . (4.17) 



Gluon-Gluon-Squark-Squark vertices 

The corresponding diagram is shown in Fig.^ (b) with its expression as follows 
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Vert (G- Gl Q{ sl , Q j m>s2 ) = ^C g, u . (4.18) 

where 



C 



1 



Hm 



48 
1 — g5s 
2 

raft 



1 + g5s 



(3 ^ 5 a6 + 8 C F {T a , T b }^ SRl(l, l) sl , s2 



(2 C F {T a , T 6 }^ (5 SRl(l, l) sl , s2 + SR2(l, l) sl , s2 ) 



+ 5*3 5 ab (7 SRl(l, l) sl , s2 + 4 SR2(l, l) sl , s2 )) 

j2 _ 9 > 



(4.19) 



Squark-Squark-Squark-Squark vertices 

The corresponding diagram is shown in Fig.0] (c) with its expression as follows 



Vert (gf M1 , Q| s2) gg )S3 , = % (Ci a* 1 ' * * + C 2 ,V ; "^ ;;! ) (4.20) 



7T 



where 



Ci = -38^2 i K1 i(- ( 42 + 4 ^ 5s ) + (78 + 8 #5 S ) N c ) S llM 5 l2jl3 

+ (6N 2 C - (42 + 4 fl 5a)) <y J1>H 5 i3 , I4 ] + ^2 (- (42 + 4 ^5s) iV c 2 - (30 + 4 gSs)) 
8n,i 2 8 im + K3 (6iV c 3 + (66 + 8 ghs) N c ) <$Z1,Z4 $12,13 
+ l_^5s^ 4 ^_ 8JV 2_ 16 ^ ^ 2 ^ 3jW +(-8iV c 3 + 32iV c ) <J IliJ4 <y,2 lB ] 

+ (12X HV - 1) [(^c + 2) 5 sltS2 5 s3M 6 lhl2 S im + (iV 3 - 4N C ) 5 slM 5 s2iS3 S llM 5 l2ii3 ] } , 
C 2 = C x (12 IA, s2 <R> sA) . (4.21) 

For brevity, we also introduce following notations 



2 

= ^2 (- R i*l,( S l,fc)- R «2,(s2,fc)-Ri*3 i ( s 3 ) fc)-R/4,(s4,fe)) , 
k=l 

K2 = ^a,(sl,l)^2,(s2,l)-R/3,( s 3,2)^4,(s4,2) + ^n,(sl,2)^Z2,(s2,2)^3,(s3,l)-^4,(s4,l), 
K3 = Rix^ s x tl )Rl2,{s2,2)Rm^ 2 )Rli,{sA,l) + ^1,(51,2)^2,^2,1)^3,(53,1)^4,^4,2), 
KA = i2n ) ( s i j i)^2,( S 2,2)-R*3,(s3,l)- R W,( S 4,2) + ^l,(sl,2)- R «2,(s2,l)^*3,( S 3,2)- R i4,(s4,l)-(4.22) 
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4.2 Effective Vertices in Mixed MSSM QCD 

In this section, we give the complete set of the non-vanishing R effective vertices in 
mixed MSSM QCD, with all internal legs SUSY QCD particles and parts of external 
legs MSSM particles. 

4.2.1 Mixed MSSM QCD effective vertices with 3 external legs 

All possible non-vanishing 3-point vertices in mixed MSSM QCD are shown in Fig.[5]. 

Scalar-Gluon-Gluon vertices 

The generic effective vertex is shown in Fig.|| (a) with its expression as follows 



Vert G; , G») = -g B* g „ ^— (± - ± £ 



]tt 2 2m w s w ysp^ op ;) 



m D I , 



v ert (h°, G; , at) = «- g ^-L-\^l + c ± £ - 1 



rn 



d » 



Vert [A or 0°, G£, G*) = 0. (4.23) 

Vector-Gluon-Gluon vertices 

The generic effective vertex is shown in Fig.[5] (b) with its expression as follows 

Vert (Z^,G",G b p ) = 5 ab £^^1-^2) — T^^q, 

Vert(A MJ G^,Gj) =0. (4.24) 



Neutralino-Gluon-Gluino vertices 

The generic effective vertex is shown in Fig.[5] (c) with its expression as follows 

Vert (Gl G\ x?) = ^ (G_fT + G+ft+) 7/i , (4.25) 
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where 



m w 



SnCpspRNt^Y, (RlR m ) ltl -4j2( RU ™ RU ™) 2 , 2 



I92m w c w s w cpsp 
+2 {RDlRD m ) 22 ) + 3m w c w c s^RN* t2 (J2 {RU^RU^^ 

m m 

-^(RDlRDrJJ + GcvCpRN^mu™ (RU} ri RU m ) 12 

m 

+6c w spRNt 3 J2™ Dm (RDlRD m ) lfi 

m 

-C* . 



(4.26) 



Scalar-Gluino-Gluino vertices 



The generic effective vertex is shown in Fig.|5] (d) with its generic expression as follows 



Vert (s, G a , = ^ (C S Q- + C^+) 



(4.27) 



The actual values of S, C s and Cf are 



C 



h° 



cf = (C 



h°\* 



wJ2 mu > (ru}ru^ -s aS pj2 m Dt { RD i RD i) 2l ) > 



a 
c 



± 



95s s^; (I>< +S5>« H«a) h ) , 



32mvi/s 
(Cf)* 



(4.28) 



Vector-Gluino-Gluino vertices 
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The generic effective vertex is shown in Fig.|5| (e) with its generic expression as follows 



Vert {y„ G a , G h ) = ^ (ClST + C v + £l + ) 7fi . (4.29) 
The actual values of V, CY and CY are 



Ct = — 



c z _ 



6 ab 

96 



2 ( (RUiRVi) i " ( RU i RUl ) 
i ^ ' 

( (rd}rd^ ^ ^ - (rd}rd, 



2.2 



-c 



8' 



ab 



E 



i 



1,1 



4sl- 1 -^3){RU]RU, 



2.2 



2s2 _l+^5£ 



3 J (rd\RD^ 



3 ) ( RDjRDi 



2.2 



cf = -ct 



(4.30) 



Neutralino-Quark-Squark vertices 

The generic effective vertices are shown in FigJ^ (/, g) with their generic expressions 
as follows 

Vert (x° n ,QhQl,s) = ^ $im {C Q 'Q- + C?fi+) , 

Vert (xl Qt, QL,s) = ^ S lm (c^CT + cf . (4.31) 

The actual values of x°, Q, Q, C_ and are 
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c: 



Ui 



a 



Ui 



c2 i 



C F [mws w spRN*iRUi t ( S)1 ) 



9G\/2mwc w s w s p 

+3m w c w s / 3RN* )2 RUi j ( Sj i) + g5s 6m Ul c w RN* 4 (RU}RUA RU l , c s+gs* ^ 

V / 1.1 ' ' 2 

—ghs &m w s w spRNn A ^RUjRUij RU^ ^+gss ^ 
+3c w m Ul RN* A RUi t{St2 )] , 

■= C F [-Am w s w spRN nA RUi t{s>2 ) 

96 y2mwc w s w s p 

+3m Ul c w RN nA RUi t ( St i) + g5s 2m w s w spRN n) i (ru} RUi^ RU^ s-g&s ^ 
+g5s 6m w c w Si3RN nt 2 ^RU}RU?j RU^ s-gs^ 

+g5s 6m UlCw RN nA (rU* RU^ RU lt , 

■= C F [m w s w Ci3RN* x-RA^s.i) 

96v 2mwc w s w cp 

-3m w c w cpRN^ 2 RDi^ sA ) + ghs 6m Dl c w RN* 3 (rd]RD^ RD^s+gs 

+g5s imwSwCpRN*^ ^RDjRDij RD^^ z+gSs ^ 

+3m Dl c w RN* 3 RD lt{St 2)] , 

—j= C F [2m w s w cpRN nA RD l>{s>2 ) 

96v 2m w c w s w c l3 

+3m Dl c w RN n>3 RDi > ( s>1 ) + g5s 2m w s w c p RN nA \ RD}RD^ RD 
—g5s Qm w c w cpRN n 2 ^RD^RDij RD t ^ 3- 8 s^ 

+g5s 6m Dl c w RN n:3 (rd}RD^^RD i 



a 



Ui 
=F 



(4.32) 



Chargino-Quark-Squark vertices 

The generic effective vertices are shown in Fig.|5| (h, i) with their generic expressions 
as follows 
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where 



Vert [Xn,UlDi 

Vert (j&Mfik) 
Vert (x-,Di,uQ 
Vert {x^Uif)Q 



Ws 



ig 2 s e 



71 z 







£jDi,Um. 




£jDl,Um. 




QUi,D m 





■fl 



(4.33) 



a 



Ui,D m>s 



a 



ci 



a 



QUi,D n 



a 



± 



Ui,D m>s 



Dl,U m ,s 



5 ij 

Um w s w c p ° F Vd ^ i~ 2m ^ CR: tl RD m , {S;1) 
+V2m Dm CR* n2 RD m ^ 2) 

+g5s 2V2m Dm CR* n<2 {RUjRU^^RD^^ 

+g5s 2V2m Ut CL n>2 (ru}RU^ RD^^ 
-g5s Am w sp CL n>1 {ru} RU^^RD m{s a^s ±) 

Um w S wSp ° F VU - Dl CL n,l RU m,(s,l) 

-V2m Um CL* n 2 RU m ^ 2 ) 

-gbs 2V2m Um CL^ 2 (rd\RD^ ^RU^^ 

5 ij r r 

64m w s w c, ° F Vu - Dl [-V2m Dl CR n , 2 RU m , {S;1) 

-ghs 2V2m Dl CR n , 2 [rd\RD^ RU m ^-^ a) 
+g5s 4m w cp CR n>1 ^RDjRD^j^RU^ 



= _ (c Ul,D 



~ [ C T ) 



(4.34) 



Scalar-Quark-Quark vertices 
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The generic effective vertex is shown in Fig.|5| (j) with its generic expression as follows 



8tt 2 

The actual values of S, Qi, Q m , v SQl ® m and a SQl ® m are 



Vert (S, Ql QQ = - C F (1 + \ HV ) (v SQlQ ™ + g5s a SQ ^"^ . (4.35) 



^UtUm _ tec a m uAm hPD&m _ ies a m Dl 5 lm 

2m w s w Sf3 ' 2m w s w cp 

a h° Ql Q m = 0> 

v H°UiU m _ iesgmu^im v H°DiD m _ iec a m Dl Si m 

2m w s w sp ' 2m w s w cp 

a H° Ql Q m = 0> 

v A°QiQ m = Q 

A°UiD m _ emu Am „A D,D m _ etpm Dl 5i m 



v H~UiD m 



v H+D t Ur, 



a H+D l U„ 



a ' "' = : , a 

2m w s w t[3 2m w s w 

v ^°QlQm _ Q 

a *°CWm = em uAlm afDlBm = 

2m w s w ' 2m H /s u , ' 



V D m ,Ul^ K(^)~ + "torn (*/»)) 

ZV Zl7l\yS w 



2\J 2m w s w 



V Vf»,DiT-7= { m u m (tp) 1 + m Dl (t p )) , 

~Vu m ,D l —7= K„(^) _1 - m Dl (tf>)) 

2\/ 2m w s w 



_ — n p 

v ,-u lDm = V \ {mui _ mD J 

2\J 2m w s w 



a *-u lDm = {mui + mDm , , 



-UiDm _ T /t , 2e 

; 2V2mvKSw 



^, 0m = fy mjJJ _ (m Um -m Dl ), 

2\/ 2m w s w 



a* DlUm = -Vv^jh—- (m Um +m Dl ). (4.36) 

2y/2m w s w 



Vector-Quark-Quark vertices 

The generic effective vertex is shown in Fig.^j (k) with its generic expression as follows 
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2 

Vert (V^, Qi, Ql) = --^- 2 C F 5 ij (1 + \ HV ) (v VQ ^ lfl + ghs a™^^.?,!) 
The actual values of V, Qi, Q m , v VQlQm and a VQl ® m are 



v AQ lQm = - ie Q Qi S lm , 
a AQtQ m = g 5 



a ZQ t Q m _ ___^Qi_fi 



iel; 

Vim 



v W-UiD m = _y\ 



it 



Dm > Ul 2V2. 



a w-u,D m 



s, 

le 



yt 

Dm ' Ul 2V2sJ 



2V2 



S, 



a w +Dl u m = V Um>Dl -^-. (4.38) 
2V2s w 



Scalar- S quark- S quark vertices 

The generic effective vertex is shown in Fig.[| (I) with its generic expression as follows 

Vert (s, Ql )Sl! Qi hs2 ) = ^ C s ^^. (4.39) 
The actual values of S, Qi s i, Qm s2 and C SQl > 3lQm > s2 are 



C h°U l>sl U m>32 = ^ Cp gij 

l92m w s w c w Sf3 
2 

v IY/io , m_c.\_ 2 , / ^„ _2 _ t 

W/3 



^ [((18 + 12 g5s) c w c a m 2 Ui + (-l) (fc+1) 4m W /m z s^s /3 s c 

k=l 

■S k A 3m w m z Si3S a+l 3) RU^ sl ^RUi^ s2 ,k)] 



Sm^Cy 



( 1 — ghs \ 

( Ac a mQ - fj,s a - c a Av t H — 8mi/,c« J RU* i{sl>2) RUi^ s2 ,i) 



( 1 — (75s \ 
+ ( 4c a mg - /x*s Q - c^c/, H - — Sm^Ca j RU^^RU^j) 
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(jA Ul }3l Um, S 2 



'Ini 



192mws w c w Cfs 



' 2 

^2 [(( 18 + 12 g^s) c w s a m 2 Di + (-l) (fe+1) 2m w m z s 2 w CpS a+/3 

k k=l 

- 4,1 3rn w m z Ci3S a+ i3) i?^* (sl fc) i?A,(s2,fe)] 



-3mj),c 



4s a m (5 - jjc a - s a A* Di H - 8m Dl s a j RD* l{sl2) RD K{s2A) 



+ ( 4s Q m (5 - fj,*c a - s a A Dl + 



A H ^ 8m A s a J RDl {sl l) RDi j{s2 , 2 ) 



= -SimTr. C F 5 lJ mu^tp)- 1 

64m w s w 

-#5s ^ + //(fy) + 8m Vl + A* v \ RU^ sh2) RU lj{s2A) 

-#5s 4m 6 + //(^) + Bmu, + A Vl \ RUl [sll) RU l:{s2 , 2) 

Q4m w s w 

-ghs 4m d + /ifo)" 1 + 1 ~ 2 95S 8m Dl + A* D }j RDf^RD^^ 



- i^-g5s Am d + //(^) 1 + 1 ^ S 8m Dl + A Dl j RD^RD^^ 

t 

64m w s w 

ghs Am d + /i(t/j) _1 - 1 - g5S 8 mUl - A* Ut ) RUf ( ^ -»RU, 



Q<t>°Ui i3l U m ,s2 _ £ 



g5s Am^ + [J,*(tp) 1 -^—8m Ul - A U{ ) RU* U , , Rl r , 



l,(sl,2) UU l,(s2,l) 
7,(sl,l) RU l,(s2,2) 



(J<t>° D l,slDm,a2 _ 



>lrn 



— ^ C F 5^m Dl (t p ) 

<6Am w s w 



g5s 4m d + (ifo) 8m Dl - A* Di J RD* l {sl2) RD h{s2tl) 

ghs 4m d + /x*^) 8m Dl - A Dl J RD* l {sl l) RD lAs2 , 2) 
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QH Ul >al b m ,s2 



1 + g5s 
Htpmut 4m (5 m{ /( 



RU l,(sl,2) RD m,(s2,2) + 

+ [(3 + 2g5s) m 2 lJi + (3 + 2#5s) vr? D j\ - m 2 w t p s w ] 



RU l,{sl,2) RD m,{s2,l) 



RU l,{sl,l) RD m,(s2,l) + 

1 - 5»5s 



V*m Dm tp + t\m Dm (A Dm - 1 + g5S 4m c 



RUl is l,l) RD m,(s2,2) \ , 



C 



H+D,„,U„ 



RD l,(sl,2) RU m,(s2,2) + 

1 — g5s 



1 + g5s 2 
fJ>tpm Dl AmQm Dl tp 



2 4m ^ ( m A + + ™ G ) + m D l ^*D l ^ 

+ [(3 + 2#5s) m^ m + (3 + 2g5s) rr? D fi - m 2 w t p s 2l3 ] 

RD l,(sl,l) RU m,(s2,l) + 



RD l,(sl,2) RU m,(s2,l) 



1 + g5s A 

V mujp + m Um A Um 4m d 



^ hS Am Dl tl (m Dl + m Um + tuq) 



(j4> Ui tSl D m}3 2 



y= C F 6 ij VX 

32\/2m w s w tB 



D m ,Ui 



R Ul( s l,l) R D m ^2.2 I ■ 

1 + g5s 



4m Dm (m Dm + m Vl + m d ) tp - A^m^tp 



RU l,(sl,2) RD m,(s2,l) 



-tp [- (3 + 2g5s) m 2 Dm + (3 + 2g5s) m 2 Vi + m 2 w c w ] RUl {sl ^RD mAs2A 
-tp rn Dm n*tp + m Dm (^-A Dm + - 

1 — g5s A , . 
4m Vl (muj + m Dm + m 6 ) 



G 



RU l,{sl,l) RD m,{s2,2) \ , 



(J<t> + £>l,slUm,s2 



Hm Dl tp + 



1 + gbs 



Am G m Dl 



32\/2rriws w ti3 F Um ' D ' | 
- 1 ^ 5s 4m Um (m Um + m Dl + m d ) - A* D m Dl 
+tf, [- (3 + 2g5s) m 2 Di + (3 + 2g5s) m 2 Um + m 2 w c w ] RD* {sll) RU m>(s2 ^ 



RD l,{sl,2) R Um,(s2,l) 



+ 



-m Um ^* + m Um ( A Um - 1 + 95S 4m d ) t p 



1 — g5s , . 
+ o 4mD ' ( mD < + mi/ ™ + m c) 



R Ul t (sl,l) R D m ,(s2,2) 



(4.40) 



Vector-Squark-Squark vertices 

The generic effective vertex is shown in Fig.[5] (m) with its generic expression as 
follows 



Vertf^Ql.!,^ 



7T 2 



(jVQl ySl Qm,s2 



The actual values of V, Qi fS i, Qm,s2 and C VQl - alQm - s2 are 



(4.41) 



+ 



96 
1 — g5s 



1 + ghs 



8SRl(l, l)sX,s2 



^ 4 {SR1{1, l) sl>s2 - SR2{1, l) sl>s2 )) ( Pl - p 2 )^ 



c zu ltSl u m , s2 = _ 6im Cf 

o IOC w S w 



1 + g5s 



,(s2,fc) 



k=l 



RUl (sl;k) RU l>ism - lQs 2 w SRU2(l,l) sl , s2 ) + (p^ 

(24i2£^ (al)1) i2C7j ) ( a2)2 ) - 12RU^ sl 2) RUi^ S 2,iy) 

-(p 2 ) M (24RUt t{sli2) RU Us2!l) - 12RU^ {slA) RU Us2j2) ))] , 



I III 



9^ (24 - 3<5 M ) RDl (sl>k) RD lj{s2>k) 
fc=i 



RDl {shk) RD Us2tk) - 8s 2 w SRD2(l,l) shs2 ) + (p^ 
{2ARD* l {slX] RD u{s2a) - \2RD* l {sl 2) RD h{s2tl) ) 

(sl,2)^A,(s2,l) - 12i?D* j ( sl)1 )i2D| j ( a2 ,2)))] 



96V2a« ° F VlDm ' Ul 



RD m ,(s2,l) 



( (pi -P2)^Y1 ( fc2 Ru i{si,k)RD mAsm ) + 

V fe=i 

"(81,2)^^,(82,1) - 8RUi {sl l) RD m ^ 2 )) - (p 2 ) M 

"( s l,l)-RAn,(s2,2) - 8i2^* (sl)2) AD m ,( s 2,l)))] , 
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r W+D ltSl U m ,s2 _ ° Vrr „ 



'■' )S 9 (pi - p 2 ) M RDl {slil) RU m>(s2 ,i) 



2 

2 

fc=l 

(4i2£) | * (sl)2) i2£/ mj ( a2l i) - 8RDi {sl l) RU m ^ s2! 2)) - ($2)^ 
(4RDl {sljl) RU m , (s2>2) - 8RDl {sh2) RU mi{s2A) ))] . (4.42) 



1 — g5s 
2 



4.2.2 Mixed MSSM QCD effective vertices with 4 external legs 

All possible non-vanishing 4-point vertices in mixed MSSM QCD are shown in Fig.^|. 

Vector-Gluon-Gluon-Gluon vertices 

The generic effective vertex is shown in Fig.^j (a) with its expression as follows 



■ 3 

Vert (V^ G°, G% G°) = ^ [CX {g^g P a + g^g™ + g^) + C 2 y w ] (4.43) 

V „„j nV 



The actual values of V, C\ and C 2 are 



C 2 A = 0, 



C Z = £ !^ Tr{T a [ T 6 j T c] j (4 44) 



Vector- Vector-Gluon-Gluon vertices 



The generic effective vertex is shown in Fig.^j (b) with its expression as follows 
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■2 



Vert (V lfl , V 2u , Gp, G b a ) = 5 ab (g^g pa + g w g av + g^g up ) 

f v ViQiQm v V 2 Q m Qi + a ViQiQ ma V 2 Q m Qi j ^ (4.45) 



where all non- vanish v v< ^ l< ^ m and a v< ^ l< ^ m are shown in Eq.( |4.38D . 



Scalar- Scalar- Gluon- Gluon vert ices 

The generic effective vertex is shown in Fig.^| (c) with its expression as follows 

■ 2 

Vert (S 1 ,S 2 ,G a p ,G b v ) = ^ 5 ab [v SlQlQm v S2QmQl - a SlQlQm a S2QmQ '^AQ) 

l,m 

where all non- vanish v s ® l ® m and a ^^ ra are shown in Eq.( |4.36( ). 

Vector-Gluon-Squark-Squark vertices 

The generic effective vertex is shown in Fig.|6] (d) with its generic expression as follows 

• 3 

Vert (V„ G a v , Q[ sl , Qi jS2 ) = ^ g, v (4.47) 
The actual values of V, Q, Q and C v ^ l - Sl ^ m -" 2 are 



C A Ql , slQm , a , = ^ T a _ q §si s2 + 2 ( 5N 2 _ g) SR1 ^ 

4 (N 2 + l) (SRl(l,l) sl:S2 + SR2(l,l) shs2 ) 



96N C 
1 — ghs 



2 



C z ^^ = 8 lm 576c l SwNc n 8 (5N 2 - 7) si SRUl(l, l) sl>s2 

+ (51 - A8N 2 ) RUl (sltl) RU l;{s2il) + 12 (2N 2 C - l) s 2 J sl;S2 
1 — g5s 



2 

2 



[-2 (N 2 C + 1) (8s* - 3) SRU2(l, l) sl , s2 



+ 12 ((2 (JV? - 3) 4 + iV c 2 4,i + 4, 2 ) RUl (shk) RU l>{s2ik) ) 



k=l 
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c 



ZD ; »i D„ 



-s, 



I I II 



57Qc w s w N c 31 I 2 



:/ K 4 (5iV 2 - 7) 4 SRDl(l, l) sl , s2 



+ (51 - 48iV c 2 ) RDl {shl) RD lXs2A) + 6 {2N 2 C - l) s 2 J sl , s2 

2 

+ 12^T (((iV 2 - 3) s 2 w + iV c 2 4,i + 4, 2 ) RDl isl}k) RD l>(s2 , k) ) 



k=l 



qW-V 1>s1 D m,s2 



-^—T-V Lui{il7 - 16N!)R U tMA} RD 



m,(s2,l) 



1 — g5s 



4^((^ c 2 4,i + M RUl {sljk) RD, 



m,(s2,k) 



k=l 



,S 



+ 2 (iV 2 + 1) (RUt (slA) RD mXs2 , 2) + RUl (slt2) RD miis2>1) ))] } 
96 J SwNc T !i V ^ i ( 17 " 16N ") RD hsi,i)RU mA s2,i) 



1 — g5s 



+ 4,2) RDl (sm RU m , 



(s2,k) 



k=l 



+ 2 (N 2 C + 1) (RDl {sl>1) RU m>{s2!2) + RDl ish2) RU m , {s2tl) ))] } (4.48) 

Vector- Vector-Squark-Squark vertices 

The generic effective vertex is shown in Fig.|5| (e) with its generic expression as follows 



Vert [V lt ^V 2u ,Ql sl , Qi s2 



Ws 



7T 



,slQm,s2 n 



(4.49) 



The actual values of V 1} V 2 , Q, Q and C^ 2 ^ 1 ^ 2 are 



(jAAQ Lal Q m , a 



(jAZQi 3l Q m>s 2 



+ 



Am 4g ^ 

1 — (?5s 



4i,«2 + 1 + ^ 16 SRl(l, l) s i, s2 



(20 SRl(l,l) sl , s2 + <iSR2(l,l) sl,s2) 



S lm C F IJ2 [Vv AQlQl [v ZQtQl + (~l) k a ZQ ^ 
\k=i 



1 — g5s 
2 

1 — g5s 



4 v A QiQi v z QiQiSR2(l, l) sljS2 



ZQiQ, 
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QW-W+U liSl U m ,s2 



C 



W-W+D, „,£». 



I,sl u m,a2 



2 2 



fe=l r=l 

+ (-l)V^) Rt^R^r) (R}Ri) k 
-J2 (l9(v ZQ ' Q ' + (-l) k a ZQ ' Qi y 

k=l ^ 

4 + (-l)(^)a^«*) 2 ) i?r, (sl>fc)J R, 



+ 



SR2(l, l) s i tS2 } , 



l,(s2,k) 



\sl,l)RU m ,( s 2,l) 



2 £ 



(> + A^ - a w+ ^) ((iiDjHDj^ _ 8 ) _ 3Slm 



1 - #5s 



W-U t D g _ a W-U t D g \ IW+D g U m _ a W+D g 



(SRD2(l, m) al , s2 + RD* W) RD mAs%2) ) } , 
~96 



— C F <{ RDl isl l) RD m:(s2A) 

RU]BU a ) ltl -%)-Z5 lm 
1 ~ g5s 4^ (v w+Dl ° 9 - a w+DlD ^ (v w ~ u s Dm _ a w-u g D m ^ 



y W-U g D m _ a W-U g D m 



QAW-Ul, sl D m ,s2 _ 



(SRU2(l, m) sl , s2 + RUl (sl;2) RU m , {s2 , 2) ) } 



144 



c v w-u lDn 



^RUl(sl,l)R D m,(s2,l) H " 



( _4 - R ^t(sl,l)-R- D m,(s2,2) + %RUl(sl )2 )RD m ,{s2,l) 
+ ^R U l,{sl,2) RD m,(s2,2))] , 



i5 ij e 
144 



c V w+ Dl u„ 



l7RD^ sl ^RU m ^ S 2,i) + 



1 - g5s 
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where the explicit expressions of v v ® l ® m and a v ® l ® m are given in Eq. (|4.38| ). 



S calar- Scalar- S quark- S quark vert ices 

The generic effective vertex is shown in Fig.[| (/) with its expression as follows 
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All of the other effective vertices which don't appear above are zero. We have 
checked that if we include all of the quarks and squarks, the terms proportional 
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to totally antisymmetry tensor e^ upcj are vanishing due to the cancelation of gauge 
anomaly in MSSM QCD. 



5. Summary 

Given the data accumulated at the LHC, searching the supersymmetry signatures 
seems feasible, which has been attracted a lot of interests from both theoretical 
and experimental sides. It is important to note that ATLAS and CMS are mainly 
focused on the MSSM to interpret what they have observed (or not observed). To 
be more specifically, a special theoretical model the so-called constrained MSSM 
was often used to reduce the number of unknow parameters from dozens to four- 
and-a-half |TUS|, |106|, |107|, |108|, |109|, |110||. Hence, phenomenologically, MSSM plays a 



special role in discovering or excluding supersymmetry. As a consequence of R-parity 
conservation, a pair of gluino or squarks are produced dominantly at the LHC and 
then cascade decays into multi QCD particles. Its irreducible background are mainly 
multijets final states, and it may also be background of other interesting processes. 
Therefore, higher order corrections especially including SUSY QCD corrections are 
predominantly in investigating SUSY phenomenology. 

In the paper, we have given a complete list of Feynman rules for the rational part 
of one-loop amplitudes in the MSSM QCD. The numerators of one-loop amplitudes 
can be simplified in four-dimensions if the corresponding rational terms are added 
correctly. It makes the NLO SUSY QCD calculation in multi-particle processes 
much easier than before. With the nature of rational part R described above, the 
regularization scheme dependence is included in the expressions of R. Therefore, 
75 problems in dimensional regularization can be clarified with the investigation of 
R in supersymmetric models. To meet the needs of practical NLO calculation, the 
results are expressed in two dimensional regularization schemes and two 75 schemes. 
The usefulness of our results will be shown in analyzing SUSY phenomenology up to 
NLO SUSY QCD corrections. 
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Figure 5: All possible non-vanishing 3-point vertices in mixed MSSM QCD. 
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Figure 6: All possible non- vanishing 4-point vertices in mixed SUSY QCD (MSSM). 
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